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Human coronavirus NL63 (HCoV-NL63) has been found in children presenting with respiratory tract 
infections (RTIs). However, the epidemiology and clinical course of this newly identified virus have not been 
fully elucidated. This study investigated the epidemiology, seasonality, and clinical features of HCoV-NL63 in 
Hong Kong children. This study consisted of two cohorts of children hospitalized in a university-affiliated 
teaching hospital. In the 12-month retrospective part of the study, reverse transcription-PCR was used to 
detect HCoVs in nasopharyngeal aspirates (NPAs). Positive samples were sequenced to confirm the identity of 
the virus and to determine its phylogenetic relationship with the HCoV-NL63 strains found elsewhere. The 
second part covered a subsequent 12-month period in which patients were prospectively recruited. Altogether, 
1,981 and 1,001 NPA specimens were studied in 2005-2006 and 2006-2007, respectively. Seventy-four (2.5%) 
HCoV isolates were identified and consisted of 17 (0.6%) HCoV-NL63 isolates, 37 (1.2%) HCoV-OC43 isolates, 
14 (0.5%) HCoV-HKU1 isolates, and 6 (0.2%) HCoV-229E isolates. HCoV-NL63 infection was more common 
in 2006-2007 than 2005-2006 (1.2% and 0.3%, respectively; P = 0.006). From 2005 to 2007, the peak season for 
HCoV-NL63 infection was in September-October, which was earlier than the peak for HCoV-OC43 infections 
(December-January). HCoV-NL63-infected patients were younger and more likely to have croup, febrile 
convulsions, and acute gastroenteritis. The majority of local HCoV-NL63 isolates were phylogenetically closely 
related to those found in Belgium and The Netherlands. In conclusion, HCoV-NL63 is an important yet 


uncommon virus among our hospitalized children with acute RTIs. 


Acute respiratory tract infections (RTIs) are responsible for 
considerable morbidity and mortality in humans. A few fami- 
lies of viruses, namely, the Paramyxoviridae, Orthomyxoviridae, 
Picornaviridae, Adenoviridae, and Coronaviridae families, ac- 
count for most of the infections (19). However, despite the use 
of a wide range of sensitive assays, the etiological agents of a 
large proportion RTIs still cannot be identified. It is not un- 
common to detect cytopathic effects in viral cell cultures but 
negative results for known viruses (24). Even in the most com- 
prehensive studies, a causative agent (either viral or nonviral) 
could be identified in only 85% of the patients (14). In part, 
this might be due to the limitations of diagnostic assays, but a 
proportion of RTIs might be caused by still unknown patho- 
gens. For instance, our recent study found that about 45% of 
nasopharyngeal aspirates (NPAs) collected from children hos- 
pitalized for RTIs were virus negative (5). 

Coronaviruses (CoVs), members of the genus Coronavirus of 
the Coronaviridae family, are enveloped viruses with a large 
plus-strand RNA genome. Three serologically distinct groups 
of CoVs have been described, and viruses within each group 
are classified by their host range and genome sequence. Up to 
2003, three human CoVs (HCoVs) were known: HCoV-229E 
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and HCoV-OC43 were identified in the mid-1960s and cause 
the common cold (11, 17, 23). The third HCoV is the recently 
identified severe acute respiratory syndrome CoV that caused 
life-threatening pneumonia with significant morbidity and 
mortality worldwide (4, 15). 

Fouchier et al. identified a novel CoV from an 8-month-old 
boy with pneumonia in 2004. Unlike HCoV-229E and HCoV- 
OC43, that virus replicated efficiently in tertiary monkey kid- 
ney cells and Vero cells. The genome sequence of this CoV was 
most closely related to the genome sequences of porcine epi- 
demic diarrhea virus and HCoV-229E. This novel virus, sub- 
sequently named HCoV-NL63 (GenBank accession number 
AY518894), could be detected in 3% of 139 patients with acute 
RTIs of unknown etiology (10). In another study, van der Hoek 
et al. described HCoV-NL63 in a 7-month-old child with bron- 
chiolitis and conjunctivitis. Five of their eight HCoV-NL63- 
infected cases were younger than 1 year old (26). All of their 
positive isolates were found during the winter of 2002, giving 
rise to a detection frequency of 7%. On the other hand, none 
of their 306 samples collected in the spring and summer of 
2003 contained HCoV-NL63. Subsequently, several groups re- 
ported that HCoV-NL63 was detected in 2 to 9% of stored 
respiratory specimens from patients suffering from a wide 
spectrum of acute RTIs (2, 3, 6-8, 18). The detection rate for 
HCoV-NL63 was higher among infants and young children (3, 
6, 8). In published studies, coinfection with HCoV-NL63 and 
other respiratory viruses was found in up to 25% of cases (2, 3). 
From those studies, HCoV-NL63 was found to be associated 
with acute RTIs ranging from the common cold and bronchi- 
olitis to pneumonia. However, the epidemiology, seasonality, 
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clinical course, and outcomes, as well as the laboratory fea- 
tures, of HCoV-NL63 infection were far from clear. Further- 
more, all except one of the published studies were based on 
retrospectively collected data (2, 3, 6-8, 16, 18). 

The objectives of this study were (i) to understand the sea- 
sonality and epidemiology of HCoV-NL63 in local children 
requiring hospitalization; (ii) to characterize the genetic diver- 
sity of HCoV-NL63 circulating in Hong Kong; and (iii) to 
delineate the demographic, clinical, and laboratory parameters 
that were associated with adverse outcomes in children with 
this emerging viral infection. 


MATERIALS AND METHODS 


Study populations. All recruited subjects were children younger than 18 years 
of age hospitalized with clinical features of acute RTI. This study was divided 
into two phases. During the first phase of retrospective case finding, 40 cases 
were randomly selected each week during a 12-month period in 2005-2006 for 
testing for HCoV infection. In the second phase of this study, 20 randomly 
selected cases were prospectively recruited each week during the subsequent 
12-month period in 2006-2007. Our interim virologic testing at the end of the 
retrospective phase could not detect any HCoV-NL63 infections among 179 
patients with wheezing illnesses and infections with respiratory viruses (influenza 
A and B viruses, respiratory syncytial virus [RSV], parainfluenza viruses types 1 
to 3, or adenovirus) identified by immunofluorescence (IF). Thus, our prospec- 
tive phase recruited only patients who were negative for these common respira- 
tory viruses. The epidemiological, clinical, and laboratory features of these sub- 
jects were collected prospectively; and their NPA specimens were tested on the 
day of collection. Blinded to the HCoV-NL63 findings, the responsible pediatri- 
cians filled in the clinical data collection form. The requests for laboratory tests 
such as complete blood count and differential and C-reactive protein level were 
made when they were clinically indicated. These parameters were captured from 
the Hospital Authority Clinical Management System currently installed for use in 
our hospital. The patients’ parents gave informed written consent, and the 
Clinical Research Ethics Committee of our university approved this study. 

Collection and processing of NPA samples. The NPA samples from the sub- 
jects were sent to the laboratory in our in-house viral transport medium, which 
contains Hanks buffer, 0.5% bovine albumin, 0.1 mg/ml gentamicin, 500 [U/ml 
penicillin, 500 [U/ml streptomycin, 2.5 g/ml amphotericin B (Fungizone), and 
NaHCO; (titrated to neutral pH) These samples were stored frozen at —80°C 
(for the retrospective study) or were analyzed fresh on the day of collection (for 
the prospective study). 

Molecular detection for HCoVs. RNA was extracted from the NPA samples 
with a PureLink viral RNA/DNA minikit (Invitrogen, Carlsbad, CA). The ex- 
tracted RNA was then transcribed into cDNA by using random hexamers with 
Superscript III RNase H-negative reverse transcriptase (Invitrogen), according 
to the manufacturer’s instructions. In view of the high degree of genetic diversity 
of CoVs, this study adopted, with modifications, the laboratory protocols for a 
low-stringency pan-CoV reverse transcription-PCR (RT-PCR) assay targeting 
the polymerase gene for the detection of HCoVs in respiratory samples (18). 
This low-stringency RT-PCR identified viral sequences that were identical or 
closely related to the published HCoV-NL63 strain (10, 18, 26). Tables $1 and S2 
in the supplemental material describe the mixture of 12 pairs of forward and 
reverse primers used for the RT-PCR assay. These primers were designed from 
the sequence alignment of 14 closely related CoVs. Table S1 in the supplemental 
material shows the primer-targeting gene positions on the sequence alignment of 
the 14 different strains of CoV starting from the 5’ position and going in the 3’ 
direction. After subtraction of the duplicate forward and reverse sequences from 
Table S1 in the supplemental material, the remaining 12 sequences were used to 
represent the sequences of the 14 different CoV strains for the design of primers 
within this region. Table S2 in the supplemental material shows the exact primer 
sequences derived from Table S1 in the supplemental material. The forward 
primers are in the 5’ to 3’ direction and the reverse primers are also reported in 
the 5’ to 3’ direction and are complementary to that of the primer positions 
shown in Table S1 in the supplemental material. PCRs were carried out in a total 
volume of 25 wl containing 2 pl cDNA, 0.25 wl HotStart Tag, 10 wM of each 
deoxynucleotide triphosphate, and 0.6 yl of each primer mixture. The samples 
were denatured at 94°C for 10 min, followed by 35 cycles of 94°C for 40 s, 55°C 
for 40 s, and 72°C for 60 s and with a final extension at 72°C for 10 min. Positive 
and negative controls were included in each PCR experiment. The PCR products 
were electrophoresed in a 2% agarose gel. The expected size of the PCR prod- 
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FIG. 1. Photograph of single and nested PCRs for 1a gene (lanes 4 
to 6) and spike gene (lanes 1 to 3) of HCoV-NL63. Lane 1, PCR with 
outer primers only (698 bp); lane 2, PCR with inner primers only (663 
bp); lane 3, nested PCR (663 bp); lane 4, PCR with outer primers only 
(839 bp); lane 5, PCR with inner primers only (525 bp); lane 6, nested 
PCR (525 bp). 


ucts was 228 bp. This pan-CoV RT-PCR assay was used to screen for CoVs in 
NPA samples collected from both the retrospective and the prospective cohorts. 

Following this low-stringency RT-PCR assay, amplification products found to 
be positive for HCoVs were purified on a MicroSpin S400 HR column (Phar- 
macia Biotech, Sweden) and were subjected to direct gene sequencing with the 
forward primer mixture (see Table S2 in the supplemental material) and an ABI 
BigDye Terminator cycle sequencing kit (version 3) on an ABI 3130 sequencer 
(Applied Biosystems, Foster City, CA). The gene sequence obtained from each 
sample was compared with the publicly accessible sequences by using the nucle- 
otide blast software of the BLASTn program (http:/Avww.ncbi.nlm.nih.gov 
/BLAST) to confirm the identity. 

Genetic diversity of HCoV-NL63 isolates. Samples positive for HCoV-NL63 
were then reamplified for the partial 1a gene and partial spike gene regions by 
nested PCR. Table S3 in the supplemental material describes the primers used 
for these PCRs (18, 26). PCRs were carried out in a total volume of 25 pl 
containing 2 wl of the previously prepared cDNA, 0.25 wl HotStart Tag, 10 ».M 
of each deoxynucleotide triphosphate, and 0.6 wl of each primer mixture. The 
samples were denatured at 94°C for 5 min, followed by 35 cycles of 94°C for 40 s, 
50°C for 40 s, and 72°C for 60 s and with a final extension at 72°C for 10 min. The 
same PCR conditions were applied to both the la gene and the spike gene 
amplifications. Positive and negative controls were also included in each PCR 
experiment. The PCR products were then electrophoresed in 2% agarose gels. 
The expected PCR products were 525 bp (1a gene) and 663 bp (spike gene) (Fig. 
1). The PCR products obtained were purified on a MicroSpin S400 HR column 
(Pharmacia Biotech) and were then subjected to sequencing on an ABI 3130 
sequencer (Applied Biosystems) by using an ABI Prism BigDye Terminator cycle 
sequencing kit (version 3) (25). 

Phylogenetic analysis. Sequence alignments of the partial 1a gene and partial 
spike gene sequences of HCoV-NL63 (10, 25) were generated by use of the 
CLUSTALW program (version 1.8) (22). Phylogenetic trees were constructed by 
the neighbor-joining method (21), and rooted phylogenetic trees were generated 
by using the PAUP* program (version 4.0 beta). The robustness of the phylo- 
genetic trees was assessed by bootstrap analysis. Bootstrap values were deter- 
mined with 1,000 resamplings of the data sets (9). Bootstrap values greater than 
70% provide significant evidence of the phylogenetic grouping. The final tree was 
obtained with the FigTree program (version 1.2.2; http://tree.bio.ed.ac.uk 
/software/figtree/). 

Statistical analysis. The results are expressed as means and standard devia- 
tions (SDs) or percentages, unless stated otherwise. Continuous variables be- 
tween the different patient groups were compared by the Student ¢ test, and 
categorical variables were analyzed by the chi-square or the Fisher exact test. 
Age, duration of hospitalization, and laboratory parameters between patients 
with and without HCoV-NL63 infection were compared by the nonparametric 
Mann-Whitney U test, owing to their likely skewed data distribution. The pro- 
portion of cases (those who were infected with HCoV-NL63) and controls (those 
who were NPA negative) and the clinical features of children infected with 
HCoV-NL63 compared with the clinical features of children not infected with 
any virus were compared by the chi-square test or the Fisher exact test. Multi- 
variate logistic regression was used to analyze the patients’ diagnoses that were 
independently associated with HCoV-NL63 infection, adjusted for age and sex as 
covariates. All analyses were done two tailed by using SPSS for Windows (version 
14; SPSS Inc., Chicago, IL), with the level of significance set at 0.05. 
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TABLE 1. Clinical and laboratory features of subjects recruited in 


2005-2006" 
Result for the following 
children: 
Feature P 
NPA positive NPA negative? 
(n = 179) (n = 1,802) 
Age (yr) 2.6 + 2.8 3:5. 316 0.001 
Male (%) 57.0 56.8 0.968 
Duration of hospitalization (days) 3.5 = 8.0 3445.4 0.835 
Duration of fever (days) Jl 2 3:0 3.7+3.4 0.044 
Clinical status (% of subjects) 
Fever (temp = 38°C) on admission 5.2 13,2 0.002 
Shortness of breath 49.3 48.6 0.921 
Cough 18.6 14.7 0.181 
Minimum arterial O, saturation 4.1 10.1 0.010 
<95% 
Abnormal chest radiograph 17.0 31.9 <0.001 
Oxygen supplementation 0 5.3 <0.001 
Laboratory results 
Hemoglobin concn (g/dl) 12.3 = 12 124+1.2 0.325 
Platelet count (10°/liter) 298.0 + 154.6 307.5 + 116.1 0.481 
Total white cell count (10°/liter) 91+ 45 12261 <0.001 
Absolute neutrophil count 4.88 + 3.70 6.81 45.31  <0.001 
(10°/liter) 
Absolute lymphocyte count 3.26 = 2.51 SAG 22.33 0.756 
(10°/liter) 
Highest C-reactive protein concn 22.1 = 31.8 36.3: = 52.0 0.019 
(mg/liter) 


“ The results are expressed as means + SDs, unless stated otherwise. 
» Cases negative for all viruses tested, including influenza A and B viruses, 
RSV, parainfluenza virus types 1 to 3, and adenovirus. 


Nucleotide sequence accession numbers. The sequence data obtained in this 
study have been deposited in the GenBank database under accession nos. 
GQ856815 to GQ856831 (1a gene) and GQ856798 to GQ856814 (spike gene). 


RESULTS 


Study populations. Altogether, 1,981 NPA samples collected 
in 2005-2006 were studied. These included 179 NPA samples 
that had tested positive for seven common respiratory viruses 
by IF, as determined by tests performed as part of our routine 
investigations for RTIs. Table 1 summarizes the subjects’ clin- 
ical and laboratory features. The mean (SD) ages of the NPA- 
positive and NPA-negative cases were 2.6 (2.8) years and 3.5 
(3.6) years, respectively (P = 0.001). NPA-negative patients 
were more likely to have a fever (P = 0.002), a low minimum 
arterial O, saturation (P = 0.010), an abnormal findings on a 
chest radiograph (P < 0.001), and oxygen supplementation 
(P < 0.001). The two groups were comparable for tempera- 
ture, blood pressure, and pulse and respiratory rates; but the 
NPA-negative cases had a higher total white cell count (P < 
0.001), absolute neutrophil count (P < 0.001), and peak serum 
C-reactive protein concentration (P = 0.019). NPA-positive 
and -negative patients were matched for the need for intensive 
care and mortality (P = 1.0 for both). One thousand one cases 
were recruited during the second prospective phase. All of 
them were negative for respiratory viruses by our routine IF 
detection method. None of these patients died or required 
intensive care. 

Detection of HCoV-NL63. Altogether, HCoVs were detected 
in 74 of the 2,982 (2.5%) cases studied: 46 from the retrospec- 
tive cohort and 28 from the prospective cohort (P = 0.792). 
Table 2 summarizes the different HCoVs detected in our sub- 
jects. The sequences of all isolates from our subjects showed 
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TABLE 2. Identification of HCoVs among subjects in 
2005-2006 and 2006-2007" 


No. (%) of subjects 


HOoV serotype 2005-2006 2006-2007 - 
(n = 1,802) (n = 1,001) 
HCoV-NL63 5 (0.3) 12 (1.2) 0.006 
HCoV-OC43 25 (1.4) 12 (1.2) 0.805 
HCoV-229E 6 (0.3) 0 0.095 
HCoV-HKU1 10 (0.6) 4 (0.4) 0.781 
Total 46 (2.6) 28 (2.8) 0.792 


“Only data for cases that tested negative for common respiratory viruses, 
including influenza A and B viruses, RSV, parainfluenza viruses, and adenovirus, 
are shown. 

» P values for all serotypes except HCoV-229E were determined by the chi- 
square test; for HCoV-229E, the results were compared by the Fisher exact test. 


=93% homology with known HCoV sequences. Five (0.3%) 
and 12 (1.2%) patients infected with HCoV-NL63 were iden- 
tified in 2005-2006 and 2006-2007, respectively (P = 0.006). 
The positive detection rates for the other three HCoVs were 
comparable in those 2 years. 

Seasonality and clinical features of HCoV-NL63 infection. 
Altogether, the combined retrospective and prospective phases 
covered a consecutive 24-month period from 2005 to 2007 and 
identified 74 HCoVs for seasonality and epidemiological anal- 
yses. Figure 2 summarizes the seasonality patterns of the dif- 
ferent HCoVs in this population. The incidence of HCoV- 
OC43 infections peaked in the winters (November to January) 
of 2005 to 2007, whereas the incidence of HCoV-NL63 infec- 
tions peaked earlier in autumn (September-October), during 
this 2-year period. HCoV-HKU1 and HCoV-229E occurred at 
a much lower incidence, and no obvious peak could be ob- 
served. May, June, and July comprised the low season for all 
HCoVs. 

Due to possible recall bias, the small number (n = 5) of 
patients with HCoV-NL63 infection in our retrospective co- 
hort were not counted in our analysis for the associations 
between this virus infection and various clinical and laboratory 
conditions. Table 3 summarizes the prospective data collected 
from HCoV-NL63-infected cases and cases negative for all 
viruses recruited in 2006-2007. Patients with HCoV-NL63 in- 
fection were younger (P = 0.007) and were more likely to have 
a discharge diagnosis of croup (P < 0.001), febrile convulsion 
(P = 0.027), or acute gastroenteritis (P = 0.028). The presence 
of HCoV-NL63, however, did not show a significant associa- 
tion with chest radiographic abnormalities or any laboratory 
parameter (data not shown). On logistic regression analysis, 
HCoV-NL63 infection remained significantly associated with 
the three diagnoses listed above (P = 0.004, 0.036, and 0.042, 
respectively). As febrile convulsion occurred only in children 
younger than 6 years of age, we also repeated the analysis for 
the significant relationship with HCoV-NL63 infection in this 
age subgroup. There was still a trend toward such an associa- 
tion (20.0% versus 4.9%; P = 0.063). The association between 
HCoV-NL63 infection and croup remained significant (8.3% 
versus 0.4%; P < 0.001) among these younger children. 

The possible associations between clinical and laboratory 
variables and infections by the other three HCoVs (n = 16) 
were also analyzed. Patients infected with other HCoVs had 
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FIG. 2. Seasonality pattern for all HCoVs and individual serologic types in NPA samples collected during the 24-month study period from 2005 


to 2007. 


shorter durations of hospitalization than NPA-negative cases 
(median, 2.0 days [interquartile range, 0.3 to 2.8 days] versus 
2.0 days [interquartile range, 1.0 to 4.0 days]; P = 0.019). The 
differences in all other clinical and laboratory parameters, how- 
ever, were insignificant between these two groups (P > 0.1). 

A total of 179 patients with acute bronchiolitis, bronchitis, or 
exacerbations of asthma in our retrospective cohort who had 
tested positive for respiratory viruses by IF were subjected to 
testing for HCoV to delineate the likelihood of any coinfection 
with common respiratory viruses. However, none of these pa- 
tients were positive for HCoV-NL63. This detection rate was 


TABLE 3. Clinical and laboratory features of HCoV-NL63-infected 
and NPA-negative cases recruited during the prospective 
12-month study period in 2006-2007% 


Result for the 
following subjects: 


Feature HCoV-NL63. NPA Pp 
infected negative® 
(n = 12) (n = 989) 
Age (yr) 0.9 (0.6-1.0) 2.0 (0.9-5.0) 0.007 


Male (%) 41.7 46.3 0.751 


Duration of hospitalization (days) 2.0 (1.0-2.8) 2.0 (1.0-4.0) 0.320 
Duration of fever (days) 2.0 (1.0-3.0) 2.0 (1.0-4.0) 0.258 
Fever on admission (% of subjects) 100 90.0 0.317 
Discharge diagnoses (% of subjects) 
Upper respiratory tract infection 33.3 34.0 0.960 
Acute tonsillitis 0 21 0.614 
Croup 8.3 0.3 <0.001 
Acute bronchiolitis 8.3 5.8 0.711 
Acute bronchitis 0 0.6 0.784 
Asthma exacerbations 8.3 6.5 0.803 
Pneumonia 0 14.3 0.157 
Acute otitis media 0 0.8 0.751 
Febrile convulsion 16.7 3.9 0.027 
Acute gastroenteritis 3a3 12.2. 0.028 
Roseola infantum 0 2.8 0.557 
Viral illness 8.3 4.7 0.552 


“The results are expressed as medians (interquartile ranges), unless stated 
otherwise. 

» Analyzed by the Mann-Whitney U test for numerical variables and the 
chi-square or Fisher exact test for categorical data. 

© Cases negative for all viruses tested, including HCoV-NL63, influenza A and 
B viruses, RSV, parainfluenza viruses, and adenovirus. 


significantly lower than that for the NPA-negative subjects 
(P = 0.019 by Fisher exact test). This finding led us to focus 
only on subjects who had tested negative for the common 
respiratory viruses in the prospective phase of this study. 

Phylogenetics of HCoV-NL63 isolates. Figure 3a and b 
shows the phylogenetic trees for the partial 1a and partial spike 
gene sequences, respectively. The phylogenetic analysis of the 
la gene sequences suggested that the majority of HCoV-NL63 
isolates circulating in Hong Kong were closely related to 
strains that have been reported in Belgium (strains BE-03-1153 
and BE-03-64880) and The Netherlands (strain NL72). This 
observation was also revealed when the spike gene sequences 
were used to construct the phylogenetic tree. A small group of 
HCoV-NL63 strains (strains CU_U0348, CU_X091, CU_P208, 
and CU_U0355) circulating in Hong Kong was separated from 
the major cluster. 


DISCUSSION 


In 2004, two groups of researchers independently reported 
on the identification and characterization of HCoV-NL63 (10, 
26). This novel virus was detected in 3% and 7% of their 
subjects, respectively. Subsequently, several groups found 
HCoV-NL63 in 2 to 8.8% of stored respiratory specimens from 
patients with acute RTIs (2, 3, 6-8, 18). In contrast, two recent 
studies suggested that this virus was less common (0.4%) than 
what we thought (16, 20). Among the published studies, coin- 
fection with HCoV-NL63 and other respiratory viruses was 
found in up to 25% of cases (2, 3). Nonetheless, none of our 
179 cases in the retrospective cohort positive for common 
respiratory viruses had HCoV-NL63 coinfection. One possible 
reason is that we have selected from the retrospective cohort 
only those cases who had wheezing disorders to test for HCoV- 
NL63, in view of the possible associations between this virus 
and exacerbations of asthma (12%) and acute bronchiolitis 
(6%), as reported in a local study (16). We might have missed 
HCoV-NL63 coinfections that occurred in children with other 
diagnoses. 
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There are four major serologically distinct groups of 
HCoVs: HCoV-OC43, HCoV-NL63, HCoV-229E, and 
HCoV-HKU1. A number of published studies showed that all 
of these viruses are capable of causing a variety of human 
diseases (6, 10, 11, 17, 23, 26). We addressed this issue by 
developing a low-stringency pan-CoV RT-PCR assay for the 
detection of HCoV-NL63 in NPA samples collected from local 
children. Recent reports on pan-CoV RT-PCRs targeting the 
polymerase and replicase la genes of CoVs helped us to de- 
velop the protocol (8, 18). This low-stringency RT-PCR al- 
lowed the identification of NPA samples containing viral se- 
quences that were identical or closely related to those of all the 
known HCoVs (10, 26). The low stringency was accomplished 
by the use of a combination of nucleotide base degeneracy, a 
low concentration of magnesium, and a low annealing temper- 
ature. Our pan-CoV RT-PCR allowed us to screen for the 
presence of known HCoVs in an efficient way. Although the 
severe acute respiratory syndrome CoV is nonexistent in our 
population at present, this pan-CoV RT-PCR would also be 
able to identify that virus. On the basis of the conserved CoV 
sequences from HCoVs that infect different species, this pan- 
CoV system offers the advantage that it may be able to identify 
previously unknown CoVs. Although we did not test for the 
sensitivity of this RT-PCR assay, Moés and colleagues re- 
ported that it was capable of detecting 5 x 10* RNA copies per 
wl of NPA sample (18). 

There were two published studies of HCoV-NL63 in Hong 
Kong. Chiu et al. reported the first local study of HCoV-NL63 
in 2002-2003 that involved 587 children hospitalized for acute 
RTIs (6). HCoVs were detected in 26 (4.4%) of the children; 
and the rates of detection of HCoV-NL63, HCoV-OC43, and 
HCoV-229E were 2.6%, 1.5%, and 0.3%, respectively. HCoV- 
NL63 could present as croup, exacerbation of asthma, febrile 
seizures, and high fever. In that study, HCoV-NL63 infections 
were noted in the spring and summer months of 2002, whereas 
HCoV-OC43 infections mainly occurred in the fall and winter 
months of 2001. Lau et al. conducted a 12-month prospective 
study in 2004-2005 of 4,181 adults and children hospitalized 
for acute RTIs in two local regional hospitals (16). In contrast 
to the methods used in our study, they performed RT-PCRs 
with primers specific for HCoV-HKU1, HCoV-NL63, HCoV- 
OC43, and HCoV-229E. They found HCoVs in 87 (2.1%) 
patients, with 0.3% HCoV-HKU1, 0.4% HCoV-NL63, 1.3% 
HCoV-O0C43, and 0.1% HCoV-229E infections being de- 
tected. Febrile seizures were more common in those infected 
with HCoV-HKU1 (50%) than in those infected with HCoV- 
OC43 (14%). The incidence of both viruses peaked in the 
winter of 2004-2005, whereas HCoV-NL63 infections occurred 
mainly in early summer and autumn. In our 24-month study 
conducted from 2005 to 2007, 74 (2.5%) patients were found to 
have HCoV infections. This incidence is comparable to that of 
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the second local study described above (16). The two studies 
had similar rates of detection of individual HCoVs and also 
found a seasonality pattern of an autumn peak in HCoV-NL63 
infections (Fig. 2). On the other hand, these findings were 
inconsistent with those reported in another local study (6) and 
in studies with Caucasians (18, 26). Selection bias is the likely 
reason that accounted for the high rate of detection of HCoV- 
NL63 in the latter studies (18, 26). Despite this, we found that 
HCoV-NL63 is associated with croup, febrile convulsions, and 
acute gastroenteritis in local children. As these diagnoses were 
not associated with HCoVs other than HCoV-NL63, our find- 
ings also suggest that the three disease associations described 
above are specific to HCoV-NL63 but not the other three 
HCoVs. Among these diagnoses, HCoV-NL63 also remained 
associated with croup when the data were adjusted by the use 
of the Bonferroni correction for multiple statistical compari- 
sons. In contrast to the other HCoVs, there was a significant 
difference in the rates of detection of HCoV-NL63 during the 
two 12-month periods (Table 2). This finding is the beneficial 
result of our two-stage study design. The many published stud- 
ies that have looked for HCoV-NL63 only in any 1-year period 
might produce inaccurate results concerning its epidemiology 
(6-8, 20). 

CoV, being an RNA virus, has a high degree of genetic 
diversity. The current literature on the sequence variability of 
HCoV-NL63 isolates in Asian countries is limited. Phyloge- 
netic analyses of the 1a gene sequences revealed that the ma- 
jority of circulating HCoV-NL63 isolates are closely related to 
strains reported in The Netherlands (strain NL72) and Bel- 
gium (strains BE03_1153 and BE03_ 64880) (Fig. 3a). Whereas 
the former strain was known from a local study (6), the strains 
from Belgium were new in Hong Kong. Besides, a small pro- 
portion of HCoV-NL63 strains (strains CU_P208, CU_X091, 
CU_U0348, CU_U0355, CU_616, and CU_743) found in 
Hong Kong were phylogenetically separated from the major 
group. For the spike gene sequences, most local isolates were 
closely related to either the prototype strain found in The 
Netherlands (strain NL496) or the strains listed above, as was 
found by la gene sequence analysis (Fig. 3b). The possible 
implications of these findings on the design of PCR primers 
and other detection assays deserve further study. We suggest 
that representative samples from both clusters of isolates be 
included when diagnostic assays are evaluated. 

The major limitation of this study is that a number of emerg- 
ing respiratory viruses, such as metapneumovirus, bocavirus, 
and rhinovirus (1, 4, 12, 13), were not detected. Apart from 
HCoVs, this study attempted to identify only common respi- 
ratory viruses by IF, and no attempt was made to identify the 
three viruses listed above. We do not know the proportion of 
HCoV-infected patients who were also infected with those 
viruses. Therefore, we cannot confidently attribute our disease 


FIG. 3. Phylogenetic tree of the HCoV-NL63 partial 1a gene (a) and spike gene (b) sequences. Bootstrap values, indicated at the nodes, were 
obtained from 1,000 bootstrap replicates and are reported as percentages. The scale bar corresponds to 0.001 substitution per site. The accession 
numbers for the reference strains are NC_005831, AY567488 to AY567494, AY675541 to AY675553, AY746451 to AY746458, AY758283 to 
AY758287, AY758297 to AY758301, DQ093116 to DQ093123, abd DQ106888 to DQ106901. CAN, Canada; Q, Australia; NL, The Netherlands; 
BE, Belgium; KR, South Korea. NL496 is the prototype sequence of HCoV-NL63. All sequences obtained from this study are named with a “CU” 
prefix. The GenBank accession numbers for our sequences are GQ856815 to GQ856831 for the la gene GQ856798 to GQ856814 for the spike 


gene. 
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associations with HCoV-NL63 to that virus alone or whether 
the diseases were caused by coinfection with HCoV-NL63 and 
other emerging viruses. In addition, the detection of a rela- 
tively small number of HCoV-NL63-infected cases in this study 
did not allow us to draw any definite conclusion on the clinical 
course of HCoV-NL63 infection and its association with clin- 
ical diseases. 

In conclusion, HCoV-NL63 is uncommon among our hos- 
pitalized children with acute RTIs, being detected in 0.3% and 
1.2% of these children in 2005-2006 and 2006-2007, respec- 
tively. The overall rate of detection for HCoVs during this 
period was 2.5%. HCoV-NL63 infections showed a peak in 
autumns (September-October) of 2005-2007, whereas infec- 
tions with HCoV-OC43, the commonest HCoV, were found in 
the winter (November to January). Compared with acute RTI 
cases in which no virus was detected, patients with HCoV- 
NL63 infections are younger and are more likely to be diag- 
nosed with croup, febrile convulsions, and acute gastroenteri- 
tis. This study suggests that the majority of the HCoV-NL63 
isolates circulating in Hong Kong are phylogenetically closely 
related to the prototype strain found in The Netherlands 
(strain NL496). Nevertheless, a small proportion of the local 
HCoV-NL63 strains were phylogenetically distinct from the 
major cluster. 
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